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• Abstract 

We apply methods of dynamical systems to study the behaviour 
^ I of the Randall-Sundrum models. We determine evolutionary paths 

' for all possible initial conditions in a 2-dimensional phase space and 

we investigate the set of accelerated models. The simplicity of our 
formulation in comparison to some earlier studies is expressed in the 
following: our dynamical system is a 2-dimensional Hamiltonian sys- 
tem, and what is more advantageous, it is free from the degeneracy 
of critical points so that the system is structurally stable. The phase 
plane analysis of Randall-Sundrum models with isotropic Friedmann 
geometry clearly shows that qualitatively we deal with the same types 
of evolution as in general relativity, although quantitatively there are 
important differences. 
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1 Introduction 



The rapid development of particle-physics-motivated cosmology represented 
mainly by superstring cosmology [1, 2, 3] has led to a change of views onto 
the standard problems of inflationary cosmology such as the (past) horizon 
and flatness problems [4, 5]. On the other hand, the astronomical observa- 
tions of supernovae la [6] strongly suggest that the universe not only had 
possibly accelerated at its early stages of evolution but it is also accelerating 
now. This puts strong constraints onto the matter content of the universe 
because only the exotic (negative pressure) matter in standard cosmology 
can lead to acceleration. However, a phenomenological nature of such an 
exotic matter in standard models expressed in terms of the perfect fluid does 
not seem to make enough connection with particle physics and that is why 
one is looking for other, more physical, descriptions. Among them the main 
proposal is quintessence or time-dependent scalar field [7] which substitutes 
an ordinary phenomenological barotropic fluid. However, there are other 
interesting proposals to express the phenomenon of acceleration which are 
related to superstring or M-theory models. The best known are pre-big-bang 
[1, 2], brane [8, 9, 10] and ekpyrotic models [11, 12]. In this paper we try 
to study the standard cosmology problems (such as cosmic acceleration and 
past horizon problems) within the framework of the brane universes. We do 
not study, for instance, the future horizon problem following the recent dis- 
cussion inspired by the S-matrix formulation within the superstring theories 
[14]. In the paper we are interested mainly in the early stage of the evolu- 
tion of the universe although the supernovae data gives restrictions onto this 
stage and they should be taken into account. 

The idea of brane universes has been first presented by Horava and Witten 
[8] who considered strong coupling limit of heterotic Eg x Eg superstring 
theory, i.e., M-theory. This limit results in 'exotic' [15, 16] Kaluza-Klein type 
compactification ofA'" = l, D = supergravity on a /Z2 orbifold (a unit 
interval) in a similar way as compactification ofA^ = l,D = ll supergravity 
on a circle results in strongly coupled limit of type HA superstring theory. 
In Hoi'ava- Witten theory there exist two 10-dimensional branes to which all 
the gauge interactions are confined, and they are connected via the orbifold, 
with gravity propagating in all 11 dimensions. After further compactification 
of Hofava- Witten models on a Calabi-Yau manifold one gets an effective 5- 
dimensional theory which has been applied to cosmology [17, 18, 19, 20]. 

Randall and Sundrum [9, 10] developed similar to Hofava- Witten scenario 
which was mainly motivated by the hierarchy problem in particle physics 
[21, 23]. As a result, they obtained a 5-dimensional spacetime (bulk) with 
Z2 symmetry with two/one 3-brane(s) embedded in it to which all the gauge 
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interactions are confined. In one-brane scenario [10] the brane appears at 
the y = position, where y is an extra dimension coordinate and the 5- 
dimensional spacetime is an anti-de Sitter space with negative 5-dimensional 
cosmological constant. The extra dimension can be infinite due to the expo- 
nential 'warp' factor in the metric 

ds^ = exp (^-2^^ [-df + d^] + dy\ (1) 

and / gives the curvature scale of the anti-de Sitter space. In the simplest 
case the induced metric on a brane is a Minkowski metric (energy momentum 
tensor of matter vanishes). However, the requirement to allow matter energy- 
momentum tensor on the brane leads to breaking of conformal flatness in the 
bulk, and the metric (1) is no longer valid. This fact is obviously related to 
the appearance of the Weyl curvature in the bulk [24, 25, 30]. The full set 
of 5-dimensional and projected 4-dimensional equations has been presented 
in Refs. [26, 27, 28]. Global geometric properties of such brane models have 
also been studied (see Ref. [29]). Generalized bulk spacetimes different from 
those of (1) have been found in Ref. [31]. Anisotropic Kasner branes have 
been immersed in AdS bulk in Ref. [32]. 

Campos and Sopuerta [33, 34] used the dynamical system methods to the 
analysis of the Priedmann, Bianchi I and V Randall- Sundrum brane world 
type cosmological models with a non-vanishing 5-dimensional Weyl tensor. 
They considered the dynamics of this model in the form of a higher-than-two- 
dimensional dynamical system. Exact analytic brane configurations with a 
vanishing Weyl tensor with perfect and viscous fiuid have been presented in 
Refs. [35, 36]. Coley [37, 38] also studied the dynamics of these Randall- 
Sundrum models - he made a step towards Mixmaster (Bianchi IX) dynamics 
and found it was not chaotic provided the matter had positive pressure. In 
this paper we show that the dynamical system which describes the evolution 
of the brane models (both isotropic (FRW) and anisotropic (Bianchi I or 
Bianchi V) types) can be represented in the simplest way in the form of 
a two-dimensional Hamiltonian dynamical system. Such visualization has 
a great advantage because we avoid the problem of degeneracy of critical 
points which appeared in the higher dimensional phase space (see Ref. [13] in 
[34]). It is well known that the existence of such critical points is a possible 
reason for the structural instability of a model. On the other hand, the 
representation of dynamics as a one-dimensional Hamiltonian flow allows 
to make the classiflcation of possible evolution paths in the configuration 
space which is complementary to phase diagrams. It also makes simpler 
to discuss the physical content of the model. Finally, the construction of 
the Hamiltonian allows to study quantum cosmology on the brane as it was 
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attempted in Refs. [39, 40] , in full analogy to what is usually done in general 
relativity [41]. 

In this paper we demonstrate the effectiveness of representation of dy- 
namics as a one-dimensional Hamiltonian flow. In this representation the 
phase diagrams in a two-dimensional phase space allow to analyze the accel- 
eration and (past) horizon problems in a clear way. We reduce the dynamics 
to a two-dimensional phase space with an autonomous system of equations. 

We deal with the full global dynamics of brane models whose asymptotic 
states are represented by critical points of the system. From theoretical point 
of view it is important how large the class of accelerated models is. We will 
call this class of accelerated models typical/generic, if the domain of acceler- 
ation in the phase space is a non-zero measure. On the other hand, if only 
the non-generic (zero-measure) trajectories are represented by accelerated 
models, then the mechanism which drives these trajectories should be called 
ineffective. Such a point of view is a consequence of the fact that, if the 
acceleration is an attribute of a trajectory which starts with a given initial 
conditions, it should also be an attribute of a trajectory which starts with a 
neighbouring initial conditions. 

Our analysis of the brane-world type model may be considered as com- 
plementary to Campos and Sopuerta's analysis [33, 34] (see also [37, 38]). 
However, it is more advantageous in many points. It is because our dynami- 
cal system is a 2-dimensional system which is the smallest possible dimension 
to study isotropic cosmological systems of equations. This allows to avoid 
huge redundancy of degenerated critical points and trajectories which appear 
in a higher dimensional phase space (cf. phase diagrams of [33, 34]) so to 
the structural instability. Finally, for a two-dimensional Hamiltonian system 
one is able to study all its properties in a configuration space rather than in 
a phase space. 

The paper is organized as follows. In Section 2 we present simple Hamil- 
tonian dynamics of the brane universes. In Section 3 we discuss cosmological 
models in configuration space, stability and the domains of the phase plane 
for which cosmic acceleration and horizon problem avoidance appear. In 
Section 4 we present phase portraits for FRW models with dust (7 = 1) 
and domain-wall-like matter (7 = 1/3) [6, 50, 51] which induces acceler- 
ated expansion (quintessence) in standard general relativistic cosmology. In 
section 5 we present simple form of dynamical systems in which constant 
coefficients play the role of the observational dimensionless density of matter 
Q parameters. In Section 6 we discuss our results. 
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2 Simple dynamics of brane universes 



The 5-dimensional Einstein equations for the Randall- Sundrum (RS) model 
with a brane aX y = location are [26, 27, 28] 

Gifj = 4)[-A(5)^J + S{y){-Xh'^:> + T(f)] (2) 

~ (51 (5) 

where Gliu is a 5-dimensional Einstein tensor, giiJ is a 5-dimensional metric. 



A(5) is a 5-dimensional cosmological constant, and h^^} — gl^J — n^n,^ is a 

4-dimensional induced metric, a unit normal vector to the brane, Tji^ 
a 4-dimensional energy-momentum tensor, A is a brane tension, and 

Stt 



,(5) 
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where the 5-dimensional Planck mass ^^^Mp is always much less than the 4- 
dimensional ^^^Mp as measured on the brane ^^^Mp -C ^^^Mp = 1.2 • lO-'^^GeV. 
This fact allows to solve the hierarchy problem in particle physics and pos- 
sibly reach the electroweak scale ~ TeV for gravity in accelerators. 
The induced 4-dimensional Einstein equations on the brane are 



^(4) 

^ 1X1/ 



where 
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(6) 
(7) 



where is the 4- velocity of an observer on the brane, A(-4') is the 4-dimensional 
cosmological constant on the brane [24, 25] and E'^,^ is the correction which 
appears from the Weyl tensor in the bulk which reads as [30, 37, 38] 



6U 



(9) 



Here lA is an effective nonlocal energy density on the brane which arises from 
the gravitational field in the bulk which is not necessarily positive and reads 



(10) 
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Since E/^i, is traceless, then its effective local pressure is p = (1/3)ZY. On the 
other hand, an effective nonlocal anisotropic stress is 

V^u - -^/«(4)A£;[^^], (11) 

while an effective energy flux on the brane is 

Q,. = -^«(4) A (E^uV + E^i^V^) . (12) 

After admission the perfect fluid with barotropic equation of state p = 
(7 — l)p, 7 e [0, 2], where p the pressure and g the energy density from (7) 
one has 

n,. = ^^V^K + ^^'(27 - i)hl^} (13) 

and the dynamics of homogeneous models in the RS brane-world scenario 
can be described by the following set of equations 

^^-^ 6 — p — 6A — ^"r+^"A4;^' 

p = -3jHp, (15) 
U = -AHU, (16) 

where H = d{\n a) / dt is the Hubble function, t is cosmological time, the brane 
tension is A > (this condition allows to recover conventional gravity) and 
«;(4) = 87rG(4) is the gravitational coupling constant (c = 1). The function U 
enters as a contribution from the Weyl tensor in the bulk. The shear scalar 
= |o-"^(Ta5 vanishes for the FRW models, whereas for the anisotropic 
Bianchi type models 

^ = -6Ha^ ^ (7^ = (7o'a-« (17) 
dt ° ^ ^ 

The first integral of the system (14)-(16) is the Priedmann equation [34] 
which reads as 



^ -—P^-^P "^ + ¥ + ^ + A4;' ^^^^ 



where k e {0,±1} - the curvature index. One can easily see from (18) that 
in the limit A — > oo one recovers general relativity. The main difference from 
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the standard general relativistic Priedmann equation is the appearance of the 
correction. This term comes as a contribution from the brane. However, 
in the general case within the framework of AdS/CFT correspondence [46] 
this correction may come from higher-derivative terms in the curvature and 
higher-dimension operators connecting the on-branc matter directly to the 
CFT [47]. 

There are two invariant submanifolds: U = and a = 0. The system on 
the submanifold U = corresponds to dynamics with vanishing electromag- 
netic part of the Weyl tensor, whereas the system on the submanifold a = 
corresponds to the case without anisotropy (FRW models). The intersection 
of these submanifolds is also an invariant. 

It is well-known that the first integral of the FRW equation can be used 
to construct a Hamiltonian function. We take advantage of this feature in 
the considered model. 

The integration of (15)-(16) gives 

p = poO~^^, U = Uoa-\ (19) 

Therefore the right-hand side of the Rayuchaudhuri equation (14) can be 
expressed in terms of the scale factor a{t) as 

(37-2)«:2 po (37-l)4)P^ 2 
— o — — CL — 7:<7oO 



6 6A 3 

a. (20) 



, A(4) 2Uo _4 
6 ■^'^(4) 



The equation (20) can be rewritten in the form analogous to the Newton 
equation of motion in the 1-dimensional configuration space {a: a E M+} 



dV , , 



where the potential function 



12A 6 6 Ak^4^ 

(22) 



and Vo = const. The plot of (22) for 7 = 1 is given in Fig.l. The equation (21) 
has a simple interpretation: the Randall-Sundrum models accelerate wher- 
ever the potential V{a) is a decreasing function of a while they decelerate 
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whenever the potential V{a) is an increasing function of a. Another impor- 
tant conclusion from (22) is that for 7 > l(p > 0) contribution dominates 
over shear anisotropy near singularity so the singularity in Bianchi models 
(especially in Bianchi IX) is isotropic and that there is no chaos [38]. In gen- 
eral relativity it is possible only for a meaningless fluid with 7 > 2 (p > p). 
The first integral of (21) is 

V{a) + ^^Vo-^. (23) 
Now we construct the Hamiltonian function 

H^^ + V{a), (24) 

and then trajectories of the system lie on the energy level H = E — const. 
The advantage of having the Hamiltonian function given by Eq. (24) is that 
one is then easily able to canonically quantize the system and allow quan- 
tum cosmological framework on the brane [39, 40] in full analogy to general 
relativity [41]. Finally, we obtain the dynamics reduced to the Hamiltonian 
flow in the 1-dimensional conflguration space 

n^^ + V(a)^0, (25) 

with V{a) given by (22) and Vq = k/2. 

Now the physical trajectories lie on the zero- energy level, 7i = E = 0, 
which coincides with the form of the flrst integral. 

The formalism can simply be generalized to the case of a general form 
of the equation of state. We use this fact and consider the equation of 
state for the (non-interacting) mixture of dust-like matter {p = 0) and an 
unknown component, labelled X, with negative pressure. It is assumed that 
dominant X-component is a perfect fluid with the equation of state specifled 
by Px — {wx — l)px and < wx < 2/3, which enables that component 
to induce accelerated expansion in standard general relativistic cosmology 
[50, 51]. The positive cosmological constant A(4) corresponds to wx = 0. 
The supernovae observations indicate that the dimensionless energy density 
of X-component Qxo ~ 0.7 is large compared to dust-like matter (pxo ^ 
— natural order of magnitude of the vacuum energy) [6] . The best-flt 
component is with wx ~ 1/3 (see paper 4 in [6]). One often represents 
the X-component in terms of a scalar field $ (e.g., inflaton) with a suitable 
potential (quintessence) [52, 53]. Conservation of energy gives p$ oc a^, 
m — — 3w$, = (w$ — l)p$, w# = const. 
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0.5 1 1.5 2 2.5 3 3.5 4 

a 

Figure 1: The potential function V{a) given by (22) for some trajectories 
depending on the initial conditions in the phase plane for k — 0, — 1,+1 
and 7 = 1 (dust). The horizontal line set by loci formed from maxima 
of all potential functions separates the regions with deceleration (left) and 
acceleration (right) (cf. Eq. (21)). Only the region with V{a) < has a 
physical sense. 
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k=+l 

k=0 
— - k=-l 



1 - 




12 3 4 



a 

Figure 2: The potential function V{a) given by (22) for some trajectories 
depending on the initial conditions on the phase plane for A; = 0, —1, +1 and 
7=1/3 (domain walls) which induces accelerated expansion in standard gen- 
eral relativistic cosmology. Similarly as in Fig. 1 there is the horizontal line of 
potential maximum points which separates to deceleration and acceleration 
regions (cf. Eq. (21)). However it appears for lower values of a. 
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The most general case of the equation of state is 

P = Px + 0= bia) -l]p= p^^3(!x-i) + i ^' ^26) 

pxo 

with wx — const, and zero stands for the pressure of the dust. 
Therefore the total energy p — Pm + Px changes as 

p = poa-' exp (- r ^^da) . (27) 



After substitution (26) into (27) we obtain 

= Pmoa-=' + pxoa-'"'^. (28) 



V{a) 



Pxo 

Then, the potential function is given by the integral 

(37(a) -2)«;2 (37(a) -l)fi;2 2 2 A(4) 2Uo _4 



ada, 



(4) 

(29) 



where p{a) is given by (28) and 7(a) comes from (26). 

After substitution the corresponding form of the equation of state (26), 
the potential V{a) can be obtained by integration of (29) 

Via) = -4^p,oa— - 4^p„oa- 
6 6 

+ I PXOPmOa — -^Pxo^ ~ 2pra0^ I 



(To _4 A(4) 2 Uo _2 k 



It is interesting that in this potential there is a term of the type pxoPmoO^^"'^"^ 
which is related to the co-existence/interaction of both matter and an un- 
known form of dark energy in the universe. 

Obviously the simplest possibility is the pure case of the unknown matter 
X such that 

P = Pxoa~^"'^, P ^ (wx - 'i-)p, 7(a) = = const, 
which gives the potential function 

V(n)- '^'^^O «-3«'x+2 '^(4) 2 ^-6wx+2 ^^0^-4 %) 2 ^0 ^-2 , ^ 

V{a) - -—pxoa - ^Pxoa - ya ' ' + 2 

(31) 



11 



i.e., wx plays the role of 7 in (22). 

As a special example let us consider the case of px = ~\Px {wx — |)- 
Then, we have 

T// N '^(4) '^(4) 2 CTq -4 ^(4) 2 ^3 , ^ /QO^ 

^(«) = --^Pxoa - ^Pxo - " "T'^ " +2" ^^'^ 

The diagram of the above function is shown in Fig. 2. Note that the 
second term in V{a) takes the form of an additive constant like the term 
k/2 and simply scales as curvature term (in general relativity curvature term 
appears for wx — 2/3 - cosmic strings). 



3 Configuration space and stability 

In our further analysis we will explore the dynamics given by the canonical 
equations. Assuming a = x and y = a — x, we have from (25) 

dy 

■ = -— = -— (34) 
^ dx dx' 

Then, we can perform a qualitative analysis of the system of autonomous 
differential equations (33)-(34) in the phase space (a, d) = {x, y) . 

Firstly, we can observe that trajectories are integrable in quadratures. 
Namely, from the Hamiltonian constraint TC = E — we obtain the integral 

r da , , 

t-to= / , (35) 

For some specific forms of the potential function (22) we can obviously obtain 
the exact solutions. However, this is not the task of our paper since we 
concentrate on qualitative analysis only. 

It is possible to make the classification of qualitative evolution paths by 
analyzing the characteristic curve which represents the boundary equation 
in the configuration space. For this purpose we consider the equation of zero 
velocity, d = 0, which represents the boundary in the configuration space. 
Because 

= -2V(a) (36) 
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the motion of the system is hmited to the region {a : V{a) < 0}. Consider 
the boundary of the configuration space given by a condition 

dM = {a:V(a)^0}. (37) 

The condition (37) together with (22) is equivalent to the equation of the 4- 
dimensional cosmological constant on the brane A(4) expressed as a function 
of a: 

\ / \ 1 ( 2 -37+2 '^^4)P0 _67+2 2 -4 ^0 -2 , o ; 

^(4)v«) = ^ ^ -«(4)Poa ^ ^ - cr^a - jj^a + 3A; 



(38) 



The plot of A(4)(a) for different k is shown in Fig. 3 (7 = 1) and in Fig. 4 
(7 = 1/3). We can see that in the latter case (7 = 1/3 there is no static 
Einstein universe. Finally, we consider the evolution path as a level of 
A (4) = const 7^ and then we classify all evolution modulo their quantitative 
properties of dynamics (compare [51, 54, 55]). We can conclude that the 
classes of admissible models in brane scenario are qualitatively the same as 
in general relativity. 

The next advantage of representing dynamics in terms of Hamiltonian is 
the possibility to discuss the stability of critical points which is only based 
on the geometry of the potential function, namely 

— if a diagram of the potential function V{a) has maxima then they 
correspond to unstable critical points; on the other hand, if a diagram of the 
potential function has minima, they correspond to stable attractors. 

— if a diagram of the potential function has an inflection point at a = ao, 
then the corresponding critical point in the phase plane is a saddle point. 

In general, the stability and the character of a critical point is determined 
by the Hessian [d'^H/dx'^dy'^]. 

In our case the Hamiltonian function (24) takes the simplest form for nat- 
ural mechanical systems (i.e., with the kinetic energy quadratic in momenta, 
and the potential energy dependent on generalized coordinates only). Then, 
the only possible critical points in a flnite domain of phase space are centers 
and saddles. 

The idea of structural stability originated with Andronov and Pontryagin 
[44]. A dynamical system S is said to be structurally stable if dynamical 
systems in the space of all dynamical systems are close, in the metric sense, 
to S or are topologically equivalent to S. Instead of flnding and analyzing 
an individual solution of a model, a space of all possible solutions is investi- 
gated. A property is believed to be 'realistic' if it can be attributed to large 
subsets of models within a space of all possible solutions or if it possesses a 
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Figure 3: The 4- dimensional cosmological constant on the brane, A^la), given 
by (38), for k = 0,-1, +1 and 7 = 1 (dust). Lines A(4) = const, gives 
quahtative classification of possible path evolutions. Let us note that there 
is a similarity to the standard FRW models. The maximum for k — +1 
corresponds to a static universe. The de Sitter model starts at the initial 
singularity. The Eddington model starts from the static universe at t — > —00 
(scparatrices on the phase plane) then reaches the singularity (such that 
a{t = 0) = Gq > 0) and evolves to infinity. The solutions are represented by 
levels of constant A(4) fixed above the A(4)(a) curve. The domain under the 
characteristic curve is non-physical. 
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Figure 4: The 4-dimensional cosmological constant on the brane, A(4)(a), 
given by (38), for k = 0, —1, +1 and 7 = 1/3 (domain walls). All models with 
A(4)(a) < and k = 0,-1 oscillate. These with k = +1 are singularity free. 
The lines of constant A(4) above the A(4) (a) curve represents the qualitative 
evolution of the model. 
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certain stability, i.e., if it is shared by a slightly perturbed model. There is 
a wide opinion among specialists that realistic models should be structurally 
stable, or even stronger, that everything should possess a kind of structural 
stability. What does the structural stability mean in physics? The problem is 
in principle open in more than 2-dimcnsional case where according to Smalc 
there are large subsets of structurally unstable systems in the space of all 
dynamical systems [42] . For 2-dimensional dynamical systems, as in the con- 
sidered case, the Peixoto's theorem says that structurally stable dynamical 
systems on compact manifolds form open and dense subsets in the space of 
all dynamical systems on the plane. Therefore, it is reasonable to require the 
model of a real 2-dimensional problem to be structurally stable. 

When we consider the dynamics of brane world models, then there is 
a simple test of structural stability. Namely, if the right-hand sides of the 
dynamical systems are in polynomial form, the global phase portraits are 
structurally stable S'^ (M^ with a Poincare sphere) if and only if the number 
of critical points and limit cycles is finite, each point is hyperbolic and there 
are no trajectories connecting saddle points. In the considered case the points 
at infinity are revealed on the projective plane. Two projective maps {z.u), 
{v, w) cover a circle at infinity given hj z = (z = oo) and v = {y = oo). 

Therefore, one can conclude that the hrane-world models are structurally 
stable. It holds because the potential function V{a) is convex up and then 
there are no non-hyperbolic centres. There are also no trajectories connecting 
saddle points. All these properties can be deduced from the geometry of the 
potential function V{a). 

Structural stability is sometimes considered as a precondition of the 'real 
existence'. Having too many drastically different mathematical models - all of 
them equally well fitting the observational data (up to measurement errors) - 
seems to be fatal for the empirical method of modern science [45]. Therefore, 
any 2-dimensional structurally unstable model is not of physical importance. 

From the physical point of view it is interesting to answer the question: 
are the trajectories distributed in the phase space in such a way that criti- 
cal points are typical or exceptional? How are trajectories with interesting 
properties distributed? For example, along which trajectories the accelera- 
tion condition, a = —dV/da > 0, is satisfied? One can easily decuce this 
from the geometry of the potential function. In the phase space, the area of 
acceleration is determined by y > or by the condition that 



(37 - 2)4po 
6 



a 



(37-1)4^"' 



a 




a-^ < 0. 



(39) 
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Prom (39) we can see that in order to obtain acceleration, the positive 4- 
dimensional cosmological constant A(4) is necessary. 

It can easily be demonstrated that if a(t) — > const as a — » 0, then the 
corresponding world model has no particle horizon. Indeed, if there exists a 
constant C such that for sufficiently large e, da/dt < C, then 



The integral on the left hand-side of this formula diverges which means that 
the time rj goes to minus infinity, and that there are no causally disconnected 
regions. Putting this in terms of variables {x,y} one needs a; — > 0, y — > 
const for the horizon problem to be solved. Now, from the Hamiltonian 
constraint we obtain that as x — > 0, then V{x) — > const. Therefore, if the 
horizon problem is solved, then the curvature effects cannot be neglected in 
the vicinity of an initial singularity. 



4 Phase plane analysis of the models with 7 = 

1 and 7=1/3 



In this section we discuss the dynamics of the brane universes in a more 
detailed way. Pirstly, it is easy to verify that the general relativistic limit 
can be recovered if one takes 




in Eq. (18). The quadratic contribution of the brane matters only when 



In order to simphfy Eq. (18) one can assume that k^^-^ = 1 or ^^^M^ — l/Sn — 
(1.2 X lO^^GeV)^ which gives 




(40) 



[31] 



p > A > (100 GeV)^ 



'^5 

and finally Eq. (18) reads as 






^A(5) + ^A ) + fa' - ka-' + 



2U 



A 



a 



-4 



(41) 
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Therefore, the dynamics depends on four independent constants. Without 
losing a generahty, wc put A(5) = Uq = po — 1. Finally, the dynamics depends 

on one parameter A (branc tension) only. 

As a proof of the effectiveness of the presented method we consider only 
two limit subcases: pure dust matter and pure matter X. 

Then, the dynamical system (33)- (34) has the form 

x = y (42) 



where A(4) = + iA. We can see that the vector field in (43) is not 

smooth. There are two important cases: ctq = 1 and ctq = 0. The former 
describe the Bianchi models (Bianchi I for k = and Bianchi V for k = — 1). 
The latter describes the FRW models {k = 0, ±1). The qualitative dynamics 
of both is similar. Due to the existence of the energy integral (25), the phase 
space is separated into the two domains by the flat model trajectory (Fig. 5). 
Let us consider first the case of the FRW dynamics with dust 7 = 1. In this 
case the effects of brane and shear in the Bianchi 1 models are equivalent. An 
interesting result which shows that there exists a maximum of scalar shear 
in the Bianchi I model was obtained by Toporensky [49] . 

The system (42)- (43) can be regularized at the origin, after introducing 
the projective coordinates 



-, u=y. (44) 

X X 



Then, we obtain 



z — —uz (45) 



The system has only three critical points in this map 

+ n 
«o = ±YY' ^0 = 0, 

Uo^O, z^Zq: 2z^ + 12z^ + Xz^ - 2A(4) A = 0. 

In order to flnd that there is one critical point with 2^0 > 0, it is sufficient to 
consider the diagram of the function 

f{z) = 2z^ + 12z^ + Xz^ - 2A(4)A. 
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Figure 5: The division of phase space for system (42)-(43) on different do- 
mains with respect to the curvature index. The flat model trajectory k — 
separates the regions of the models with negative and positive curvature. 
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u 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 



Figure 6: The phase portrait of the system (42)-(43) with A = 10*^ and 
7 = 1. The Einstein-de Sitter and Eddington models are represented by 
separatrices in the phase plane. In the neighbourhood of separatrices one 
can see the Lemaitre- Eddington (L-E) evolution with characteristic quasi- 
static regime. The closer to a critical point the trajectory is, the longer the 
time of a quasi-static stage evolution. All phase curves lie on algebraic curves 
given by the first intergral (Hamiltonian constraint). The acceleration region 
is situated to the right of the saddle point, therefore for the L-E universes 
the acceleration begins in the middle of quasi-static stage. 
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Figure 7: The phase portrait of system (42)-(43) with A = 10^ and 7 = 1/3. 
The negative curvature term and the brane term in potential (22) are of the 
same type. There are two types of trajectories which expand to the maximal 
scale factor a and then recoUapse or vice versa. In the first type of evolution 
the closed models start with the singularity a(0) = and d = 00, reach the 
maximum radius and recollapsc while the other class of models recollapse 
from a(0) = 00 and d = — 00 to the minimum value of a, and then expand 
to infinity. As in the 7=1 case, the acceleration region is on the right to 
the saddle point. 
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If A(4)A > then /(O) < and because f{z) is strictly increasing as z > 
then there exits only one critical point such that zq > 0, such that f{zo) = 0. 
Let us note that for A(4) < there is no such a point. The phase portrait of 
this system is presented in Fig. 8. 

For the completeness of the analysis of the dynamical system, it is nec- 
essary to consider another projective map {v, w) 

v = w = -. (47) 

y y 

After regularization and after introducing a new time variable rj: drj — dt/w^ 
we obtain from (42)- (43) and (47) 

v = v( Iv'w' + - + ^vW] (48) 

yb OA O A J 

w^w(w' + Iv'w' + ^v' - + ^vW) . (49) 

\ D OA o A J 

In this map we can find only one critical point v = and w = 0. As it was 
mentioned before, it is interesting to observe how much time the trajectories 
spend in the area of accelerated expansion defined as 

{{x,y):x > Ay > 0} 

{X, y) : ^x-- + + la,.-" - ^ + l^-' < 

For the special case of dust we obtain in ti)-coordinates 

{1 1 2 A. 2 1 

(z, u) : -z^ + — + -ao^' - ^ + -^^>0Az>0|. 

The dynamical system for pure accelerating matter has the form analo- 
gous to (42)-(43) but now wx coincides to 7. As an illustration we consider 
the case of px = —\px (domain-wall-like matter), i.e. wx = \- Then, we 
have 

x^y (50) 

y=-- Ix-' + ^x. (51) 
^ 6 A 3 ^ ' 

In this case, the term which comes from the brane and which is proportional 
to vanishes (see Eq. 18), and only the term which comes from a non- 
vanishing Weyl tensor is present in (50)-(51). 
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X 

Figure 8: The phase portrait of system (45)-(46) with A = 10^ and 7 = 1. 
All closed generic models are of two types. The first starts from the anti- 
de Sitter and finishes at de Sitter stage. The second starts from an initial 
singularity, expands to a maximum size, and then recoUapses to a second 
singularity. Non-generic (exceptional) cases corresponds to the separatrices 
going in or out of the saddle point, which represent the static universe. For 
all open and expanding universes the de Sitter model is a global attractor, 
and for open and contracting universes the anti-de Sitter is a global repeller. 
The singularity is reached at 2; = 00 and the half-plane 2; < has no physical 
sense. All points at the infinity which arc represented by {z ~ 0} axis are 
hyperbolic, therefore the system is structurally stable. The same result is 
valid in the plane. The brane effects produce p^/A term which is 

formally equivalent to the effects of the stiff matter equation p = p (or a 
massless scalar field). 
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Figure 9: The phase portrait of system (52)-(53) with A = 10^ and 7 = 
1/3. The acceleration region is situated on the left from the saddle point. 
Therefore the recoUapsing and then expanding models lie permanently in 
the acceleration domain. The Eddington model is also in this region. The 
Lemaitre-Eddington type models start accelerating in the middle of quasi- 
static phase. 
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After introducing the projective coordinates 



1 y 

z — —, u — — 

X X 



we obtain 



z — —uz (52) 
u=lz-lz' + ^-u\ (53) 



The critical points are 



/A (4) 

^°^^V~3~' ^0 = 0, 

1 2 -A. 

Mo = 0, Zo = z: f{z) = -z--z^ + ^ = Q. 

It can be shown that /(-Zq) = always exists because /(O) > and /(oo) = 
— oo. 

In turn, in the map 

1 X 
V — -, w — — 

y y 

we have 

s^v fl 3 , A(4) 4 2 A , ^, 

dry V6 3 A y ^ ' 

dw ndw 3 /I 3 , A(4) , 2 A 

— = W^— = w'^-w\ -vW^ H —w^ - -V* . (55) 

dr) dt \6 3 A y ^ ^ 

Prom (50)- (51) one can easily determine the domain of acceleration which 
corresponds to the domain of x in which V{x) is a decreasing function of its 
argument 

_dV 
dx 

where x^ax : %\x=x^^^ = or y{x^^^) = 0. It can be checked that y is a 
strictly increasing function of a (y > 0) and y(0) = — oo, so there is a single 

point Xmax- 

From (51) we can also observe that as x goes to zero then y goes to 
infinity, which means that the model has a particle horizon (A(4) is negligible 
and then y (x x~^). It is due to the presence of 7^ 0. 
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Prom (42)-(43) and the integral of energy we can see that we have the 
horizon in a generic case if po 7^ A 7^ A ctq 7^ 0. In the special case 
of Mo = "^o = we have a model without a horizon (the FRW models with 
vanishing electromagnetic part of Wcyl tensor). 

From (50)-(51), we obtain that as a; — > 0, then y = ^, i.e. y oc t and 
x{t) oc t'^. 

In general, from the first integral (25) one observes that the generic evo- 
lution of cosmological models without a horizon is when V{a) — > const, (zero 
is also possible) as a ^ 0. 

Effects of A (4) are always negligible near a singularity and if Mq 7^ we 
have always asymptotically V{0) — 00. For uq — ao — and 7 < | (domain 
walls on the brane), then V{a) goes to a constant and there is no horizon in 
the model. Then, the generic behaviour near the singularity is x(t) oc t, i.e., 
it is Milne's evolution (in general relativity this appears for 7 = 2/3 (cosmic 
strings on the brane [51]). 

5 Simple dynamics in terms of density pa- 
rameters Q and the exact solutions 

Independent observations of supernovae type la, made by the Supernovae 
Cosmology Project and the High z Survey Team [6] , indicate that our Uni- 
verse is currently accelerating. There is a fundamental problem for theoretical 
physics to explain the origin of this acceleration. If we introduce the cosmo- 
logical constant A(4) and assume that the Universe is flat, then the best-flt 
model is for the cosmological constant density parameter f^A(4),o — 0.72 and 
for the dust density parameter equal to Q^_o — 0.28 (index "0" refers to the 
present moment of time). 

Our formalism gives natural base to express dynamical equations in terms 
of dimensionless observational density parameters Q and to compare the 
results with supernovae data. However, before we study these quantities 
in detail following the discussion of Refs. [51, 55, 56, 58] we introduce the 
notation in which it is easy to tell which models are exactly integrable. The 
Friedmann Eq. (18) with the help of the conservation equations (19) and 
(20), can be rewritten in the form 

lfday_CGR^Cx k AaCu 

'^['dij ^ + T + ^' ^^^^ 
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where we have defined the appropriate constants {k^^-^ — SttG) 

Cgr = ^a^^^, (57) 

C, = ^a^-^^ (58) 
Cu = -^a^U, (59) 

and Cgr is a of general relativistic nature. It is easy to notice that the follow- 
ing cases can be exactly integrablc in terms of elliptic functions [55, 51]: 7 = 
(cosmological constant), 7 = 1/3 (domain-wall-like matter) and 7 = 2/3 (cos- 
mic strings). The first case is the easiest since in this case the first two terms 
on the right-hand-side of (56) play the role of cosmological constants similar 
to A(4). The next two cases involve terms which were already integrated in 
the context of general relativity. For 7 = 1/3 (domain- wall-like matter on the 
branc) the general relativistic term with Cqr in (56) scales as domain-walls 
in general relativity while the term with Cx scales as cosmic string (curva- 
ture) in general relativity. For 7 = 2/3 the term with Cx scales as radiation 
in general relativity. Then, the problem of writing exact solutions reduces to 
the repetition of the discussion of Refs. [55, 51]. We will not be doing this 
here. For other values of 7 = 4/3; 1; 2 the terms of the type and 
appear and the integration involves hyperelliptic integrals. 

Coming back to observational quantities we now define the dimensionless 
observational density parameters [56, 57, 58] 

.2 



K 



^GR = (60) 

2 



where the Hubble parameter and the deceleration parameter read as 

H = (65) 
a 

q = (66) 



27 



so that the Priedmann equation (18) can be written down in the form 



^GR + + + f^A(4) +^u = l- (67) 

Using (60)-(66) the equation (20) can now be rewritten as (compare Eq.(lO) 
of [56]) 

On the other hand, it is useful to express the curvature of spatial sections by 
observational quantities using (67) and (68) 

^ ^nGR + ^^^n^ + 2nu-q-i. (69) 



H^a^ 2 2 

These relations (68) and (69) together with the exact solutions (which here 
are certainly available for 7 = 0; 1/3; 2/3) (see Refs. [58, 59]) are useful 
in writing down an explicit rcdshift-magnitude formula (generalized Hubble 
law) for the brane models to study their compatibility with astronomical 
data from supernovae. Such a comparison will be presented elsewhere [60] . 

In this paper we only study the formalism which allows to formulate 
our phase space quantities in terms of the observational parameters fl. Let 
us consider a brane universe filled with an unknown type of matter with 
barotropic equation of state Pi = (wj — l)pi. Then, it is useful to rewrite the 
dynamical equations to a new form using dimensionless quantities 

x = —, T= \Ho\t. (70) 
ao 

The basic dynamical equations are then rewritten as 

i 

x^^Yl^,,o{2-Swi)x'-'^% (72) 

i 

where = {^gr, ^x, ^u, ^Ai)- 

The above equations can be represented as a two-dimensional dynamical 
system 



y=lj2^,,o{2-3w,)x'-'^^. 



28 



As we have mentioned already the dynamics can be always reduced to the 
form of the Priedmann models of general relativity with a certain kind of a 
non-interacting multifluid. 

Of course, the dynamical system (71)-(72) has the Hamiltonian 

n = ^ + V{x), (73) 

where V{x) = —^flk,o — | ^1,02^^"^"'% which should be considered on the 
zero-energy level. 

As an example of application of these equations to study the problem of 
cosmic acceleration, consider the case of D,u,o — and f^A(4),o 7^ 0. It emerges 
that at present our Universe accelerates provided that 

Qx,o(2 - 3wx)x^-^'"'' + ^^A,o(2 - 6wx)x^-''"'' - 2Q^^^^,oX > 0, (74) 

where wx refers to any type of matter X with the equation of state px — 
{wx — ^)qx (for f^A^^j = it is possible if wx < 1/3)- 

It is convenient to introduce a new variable z — x'^'^^ > 0, then Eq. (74) 
reduces to the quadratic inequality 

fiA,o(2 - 6wx)z^ + nx,o{2 - 3wx)z - 2fiA(,,,o > 0. (75) 

In the special case of the component X in the form of the wall-like-matter 
Wx = 1/3, the Universe accelerates if 

nx,oz > -2f^A(4),o, (76) 

i.e., if A(4) > 0, Wx > 0, the universe always accelerates. 

In general, the domain of acceleration depends on the solution of the 
inequality (75). We assume that we are in the region of acceleration i^A(4),o > 
0, Wx > and consider two cases. 
Case 1: wx > 1/3. 

There are always two solutions for z with opposite signs. One solution is for 
< z < z+ (or x+ < X < 00), where x+ is given by 



{3wx - 2)Qx,o + J[{2 - 3wx)nx,o? - 16(1 - 3w;x)f^A(4),o^^A,o 

X4- ^ 



{3wx - l)^x,o 

' (77) 

In this case we have a minimum value of acceleration at 

(2 -3«;x)n^,o^ 



•^'min 



2{6wx - 2)Qa,o 



(78) 
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Case 2: wx < 1/3. 

The universe accelerates for every x. 

Let us note that in Case 1 the minimum value of x can always be expressed 
in terms of redshift z = x^^ — 1 and from the inequality < 1 we obtain 
an additional restriction on parameters fij^o- 

6 Conclusion 

In the paper we studied the dynamics of Randall-Sundrum brane universes 
with isotropic Friedmann geometry. Our approach is the simplest in that we 
formulate the problem in a 2-dimensional phase space and not in a higher- 
dimensional phase space as it has been done recently. First of all, this is 
the smallest possible dimension to study the isotropic cosmological systems 
of equations. Then, it allows to avoid the degeneracy of critical points and 
trajectories which appear in a higher dimensional phase space and so to avoid 
structural instability of these models. Also, since the system is Hamiltonian 
one is able to study its properties in a one- dimensional configuration space 
rather than in phase space. In such a space the motion of the universe point is 
reduced to the motion of a particle with the potential energy V{a), where a is 
the scale factor (see Section 2). On the other hand, Hamiltonian formulation 
can easily be applied to quantum cosmology on the brane. 

In Fig. 1 and 2 we can see that there is no qualitative difference in the 
shapes of the potential function V{a) as 7 varies for cases under studies. Only 
the values of a for which the potential is negative have the physical meaning. 
The maximum of the potential function corresponds to an unstable saddle 
point on a phase diagram (a critical point solution of the zero energy level). 
The de Sitter solution is a point at infinity. Near the initial singularity 
we have a solution a{t) oc t^^^'^ which differs from a general relativistic one 
a{t) oc ^2/37 . 

From the theory of qualitative differential equations in Sections 3 and 4 
we obtained the visualization of the system evolution in the phase plane {x. x) 
and analyzed the asymptotic states and concluded that the brane models are 
structurally stable. We studied the solution of the (past) horizon problem 
and the initial conditions for acceleration of our Universe. We have the neat 
interpretation of a domain of acceleration as a domain in configuration space 
where the potential function decreases. Therefore from the observation of 
the potential function V{a), we can see the acceleration domain a > amax, 
with y(amax) = 0, which is independent of the curvature index k. On the 
other hand, the (past) horizon problem is solved when V{a) — > const (may 
be equal to zero) as a ^ 0. If we find trajectories for which y — a goes to 
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infinity as x — a ^ 0, then the horizon is present in such a model. This is the 
case of dust brane model as well as a brane world filled with domain- wall-like 
matter with 7 = 1/3. 

In general our model does not solve the horizon problem due to the ex- 
istence of matter dominated phase of the early evolution where the term 
dominates and so H ex p oc a"^'^ {a{t) oc t'^^^'^). Only in the special case 
of heavy domain- wall- like matter (7 < 1/3) the horizon problem is solved, 
because the evolution near the singularity dominated by matter (brane ten- 
sion) is a oc a^~^'^ . However, the existence of such a singularity requires 
U = (Weyl tensor contribution from the bulk vanishes). Therefore, the 
case 7 < 1/3 is strictly distinguished because as x — > then y — > const. For 
7 = 1/3 the term p^/X vanishes and brane effects are negligible. 

In the phase portraits we can observe similarities as well as differences 
in both considered cases. In the dust-like matter case there is a quasi-static 
stage first discussed by Lemaitre (called 'loitering stage' in [48]). Such quasi- 
static stages, present in the case of wall-hke-matter, corresponds to configu- 
rations in the vicinity of a critical point. The acceleration does not depend 
on the value of Qk, but in the phase space there are different sets of initial 
conditions which provide acceleration. 

We conclude that the classes of admissible models in brane scenario are 
qualitatively the same as in general relativity, i.e., there exists a homeomor- 
phism in the phase space which transforms trajectories of brane models into 
general relativistic ones without a change of orientation. However, quanti- 
tatively the models differ from general relativistic ones - for instance, they 
asymptote the singularity in a different way. 

Finally, we reformulated the phase space quantities interms of the obser- 
vational parameters Q and presented the formalism in which one is able to 
perform the observational tests (redshift-magnitude relation) of brane models 
which we will address in a separate paper. 
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